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Abstract 

We describe the cosmological dynamics of perfect fluids within the framework of effective 
field theories. The effective action is a derivative expansion whose terms are selected by 
the symmetry requirements on the relevant long-distance degrees of freedom, which are 
identified with comoving coordinates. The perfect fluid is defined by requiring invariance 
of the action under internal volume-preserving diffeomorphisms and general covariance. At 
lowest order in derivatives, the dynamics is encoded in a single function of the entropy 
density that characterizes the properties of the fluid, such as the equation of state and 
the speed of sound. This framework allows a neat simultaneous description of fluid and 
metric perturbations. Longitudinal fluid perturbations are closely related to the adiabatic 
modes, while the transverse modes mix with vector metric perturbations as a consequence of 
vorticity conservation. This formalism features a large flexibility which can be of practical 
use for higher order perturbation theory and cosmological parameter estimation. 



1 Introduction 



The success of effective field theories for studying physical systems comes from the fact that they 
allow to capture in a single picture the universal long distance (low energy) properties of models 
that are instead intrinsically different at much shorter scales. Since cosmological problems are 
often characterized by well separated scales, the language of effective field theories provides a 
powerful tool for the study of cosmological evolution, in particular inflationary and dark energy 
dynamics. Most of the works on cosmological evolution that are based on effective theories can 
be broadly classified in two different categories, depending on whether they aim to describe the 
full (effective) action [1-6] or, conversely, focus on the dynamics of the perturbations around some 
background [7-29]. The formalism we develop in this work pertains to the former class but, as we 
will see, gives a general and straightforward effective expansion for cosmological fluid perturbations 
on a Friedmann-Lemaitre-Robertson- Walker (FLRW) metric (and may also be extended to other 
backgrounds). 

There has recently been a renewed interest in the understanding of fluid dynamics from an 
effective field theory point of view [30-36] . This general approach is based on the identification of 
the relevant long-distance degrees of freedom and their symmetries, possibly including spacetime 
symmetries such as Galileo, Poincare, or diffeomorphism invariance. The dynamics is organized 
in a systematic derivative expansion that makes the theory predictive at low energies. This 
formalism (or variations of it) has already been applied to describe some aspects of perfect fluids 
and superfluids [33-35]. However, although fluids are ubiquitous in cosmology, the cosmological 
applications of this framework have been, to the best of our knowledge, largely neglected 1 . In this 
work we describe the basic set-up for perfect fluids, adopting the general principles of effective 
field theory to describe long-distance relativistic dynamics in cosmology and their interplay with 
metric perturbations. 

Neglecting chemical potentials, we consider perfect (dissipationless) fluids that carry no con- 
served charge. We thus focus on minimal fluids that involve only three degrees of freedom asso- 
ciated with the position in space of a fluid element. These degrees of freedom can be identified 
with comoving coordinates and, as we will see, the symmetries of the fluid determine the oper- 
ators that appear in the action. In particular, perfect fluids are invariant under internal spatial 
diffeomorphisms that preserve the volume. As a result, the relativistic dynamics is fully described 
at lowest order in derivatives by a single operator. By construction, these fluids support both 
longitudinal (compressional) and transverse (vortices) excitations, each with its own dynamics. 
We discuss perturbation theory for a fluid with an arbitrary equation of state, and describe the 
coupling to scalar, vector and tensor metric perturbations. Vector metric perturbations mix with 
vortices, while adiabatic perturbations correspond to the compressional modes of the fluid. Our 
formalism can cover and extends other frameworks (e.g. P(X) theories, where the dynamics of 
fluid vortices is instead absent). 

Following the formulation developed in [30-33], which can be connected to earlier works [38-43], 
we extend the effective theory of perfect fluids to generic metric backgrounds in Section 2. We then 
focus on the FLRW case in Section 3, and move on to discuss matter and metric perturbations 

1 See however [37] for an application of this formalism to a Lorentz violating dark matter model. 
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in Sections 4 and 5. Adiabatic perturbations are introduced in Section 6; and the conclusions 
and future work directions are presented in 7. In the Appendix A we review perfect fluids in the 
context of relativistic hydrodynamics. In the Appendix B we comment on the relation between 
Eulerian and Lagrangian formulations for fluids, which we both use through the text. We discuss 
vorticity conservation using ADM variables [44] in Appendix C. 

2 The effective theory of perfect fluids 

A perfect fluid is described by a functional of three spacetime scalar functions $ a G R with 
a G {1, 2, 3} , that define, at any time, an isomorphism between the three-dimensional coordinate 
space of an observer O and a continuum of points T (the fluid) [30,39]. These functions $ a label 
a generic fluid element 2 in J 7 , whereas the map $ a — > x l (r, <E») gives the position in real space 
of the fluid element $ at a given time r. In other words, x % = x l (r, $) is the trajectory of the 
fluid element $. This corresponds to a Lagrangian description of the dynamics, see Figure 1. 
Vice versa, the inverse map x % — > $ a (r, x) returns the fluid element $ that is sitting in x at the 
time r, providing an Eulerian description of the dynamics. In a stationary background state, the 
isomorphism can be chosen such that $ a (r, x 1 ) = x a . 

Since both O and J 7 are isomorphic to 1R 3 , there is always a change of spacetime coordinates 
O — > O such that x % i— >■ x l = $ l (r, x^) so that $ a can be naturally interpreted as the comoving 
coordinates of the fluid. This means that their variation along the fluid four- velocity u^ 1 = dx^/dr] 
(being rj the proper time) is zero 

u»d^ a (x,r)=0, «% = -l. (2.1) 

As we will see next, these two conditions characterize completely the fluid four- velocity, once the 
symmetry properties of the system have been chosen. Indeed, the whole structure of the effective 
theory is fully determined by the symmetries of the fluid. 

Since the fluid must be homogeneous and isotropic, the internal coordinates have to satisfy 
the symmetries 

$« — + c « (2.2) 
$ a — >R a b §\ ReSOii) (2.3) 

where c a and the matrix of elements R% are constant in space and time. In a homogeneous 
and isotropic model of the universe, even though the background solution $ a (a; l ,t) = x a , which 
represents the ground state of the system, spontaneously breaks spatial translations and rotations, 
the diagonal combination of internal (acting on $) and space (acting on x) symmetries is left 
unbroken [30]. These diagonal symmetries ensure that the perturbations (or in other words, the 
excitations of the fluid) 7r a = $ a — x a propagate in a homogeneous and isotropic background. 

2 We are focusing on mechanical fluids that at each point in space are fully described by three degrees of 
freedom. See Appendix B for more details and the relation between these fluids and the continuum limit of 
relativistic uncharged point-like particles. 
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Figure 1: The map between <3> and x coordinates is depicted. At any given (conformal) time r , a fluid 
element labelled by <3? occupies a position given by x(t, <&). If the inverse function is considered, any 
spacetime point (r, x) is mapped to a fluid element In this picture, the <3? coordinates are scalar fields 
of spacetime. 



In addition, we demand invariance under volume preserving spatial diffeomorphisms 

^ —> , det = 1 (2.4) 

that distinguish a perfect fluid from a gel (or jelly), which is a homogeneous and isotropic solid 
[30]. It is clear that the spatial volume preserving diffeomorphisms (2.4) include the symmetry 
transformations (2.2) and (2.3), but we have highlighted the latter two because of their clear 
geometrical meaning. 

The low energy effective theory for a perfect fluid is given by a Lagrangian £ organized in 
a derivative expansion. Since the dynamics is invariant under spacetime diffeomorphisms and 
the transformations (2.2), (2.3) and (2.4), the lowest order Lagrangian must be a function of the 
determinant of the matrix B, whose elements are given by [30, 38] 

B ab = g^d^ a d v $ b . (2.5) 
The effective action that describes the low energy dynamics of the fluid is then 

S m [$] = I d A x^-gC m (2.6) 



where the Lagrangian (density) is a function of the determinant of B 

Cm =F(b) (2.7) 
b =Vdet B (2.8) 

that encodes the long distance properties of the fluid. We are implicitly assuming that there are 
no extra symmetries that could forbid any possible Lagrangian F(b). Otherwise, the low energy 
Lagrangian would start at the next order in derivatives: g^u^d^b, where g(b) is an arbitrary 
function of b [32]. However, this term can be recast into a higher derivative term by a field 
redefinition [32]. In the following we will consider only the lowest order Lagrangian (2.7). 
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The equations of motion that come from (2.7) are 

d, [V=gg» u bF b (B- 1 ) cd d v ^ d \ = (2.9) 

where F b denotes the derivative of F(b) with respect to b. The gravitational energy- momentum 
tensor of the system is 

Tfxv = - J—^Z = 9,uF - bF b {B- l Y d d^ c d^ d . (2.10) 

This corresponds to the energy-momentum tensor of a perfect fluid 

Tpv = (p + V)UpU v + pg^ (2.11) 

whose components can be easily identified. The conditions (2.1) determine the fluid four- velocity, 
which is 

' e^e abc d a ^ a d^ b d^ c . (2.12) 



6b^g~ 

Contracting T^ u with u^u v (and recalling that u 2 = — 1) one gets 

p = -F (2.13) 



and therefore the fluid is not only perfect but also isentropic [40] . For a physical interpretation of 
£ m = —p see Appendix B. Finally from the trace we also obtain the pressure 

p + p = -bF b . (2.14) 

The fluid is thus barotropic because the pressure depends only on the energy density given that p 
and p are both functions of b only. 
We also find that 

V = {B' l ) ab d^ a d^ b = g, u + Ulx u v (2.15) 

is the standard projector on hypersurfaces orthogonal to the four-velocity of the fluid. 

Notice that for dust (pressureless matter) p = and so F oc b, while for radiation p = 3p 
and F oc 6 4//3 . In general, the relation between the background energy density and pressure is, as 
usual, given by the equation of state, which is defined as 

p = wp (2-16) 

and therefore using (2.13) and (2.14) this gives 



w = -l + b = _i_|_ — 2.17 
F a log b 
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The current 3 

jn = _ bu n (2.18) 
is covariantly and identically (i.e. off-shell) conserved 

J% = (2-19) 

and it is identified with the entropy current [32,36]. The comoving entropy density J^u^ = s 
becomes 

» = 6=^ (2.20) 
where T is the fluid temperature, given by 

T — —F b . (2.21) 

The invariance under spatial diffeomorphisms (2.4) gives rise to an infinite set of (on-shell) 
conserved currents 

Jfc = -bF b (B-^d^e ^) Jfc ;M = (2.22) 
where e a ($) is an arbitrary transverse function of $ 

^($) = ^=0. (2.23) 

This constraint, generically solved by e a ($) = e abc dbf c (^), ensures that $ a — > $ a + e a ($) is an 
infinitesimal volume preserving diffeomorphism. In the Appendix C we explicitly construct the 
associated conserved charges. 



3 Effective perfect fluids in FLRW 

We have seen that the perfect fluid form of the energy-momentum tensor (2.10) comes from impos- 
ing the conditions of homogeneity, isotropy and invariance of the action under volume preserving 
diffeomorphisms of the internal coordinates. In cosmology, the first two assumptions directly lead 
to the FLRW metric. In fact, the FLRW metric is a purely geometric consequence of requiring 
that the universe appears homogeneous and isotropic to free falling observers [46,47]. For sim- 
plicity we assume from now on that the background metric of the universe is of flat FLRW type 
and work with conformal time r, 

ds 2 = -a 2 (r)(rfr 2 + S^dx*) . (3.1) 

3 We fix the normalization of such that the temperature in (2.21) matches the standard thermodynamics 
normalization, T = (dp/ds). 
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The usual equations of motion that govern the background cosmology in the metric (3.1) simply 
read: 

U 2 = —a 2 p T H = —a 2 {p T + 3p T ) (3.2) 

where H = a/a. Newton's gravitational constant is denoted by G and the quantities px and px 
represent the total background energy density and pressure if several fluids are present. Let us 
point out that if there are several fluids that only interact among themselves through gravity, each 
of them satisfies the following Friedmann equation for the background 

p<x + 3(1 + w a )Hp a = . (3.3) 

From the equations (3.2) and the expressions (2.13) and (2.14) one sees that all self-accelerating 
solutions (which require wt < —1/3) have to satisfy the condition 

3bF Tb < 2F T (3.4) 

where Ft is the function (2.7) that describes the total energy density of the fluid admixture. 
Notice also that (2.17) implies that a perfect fluid will have an equation of state that is close to 
that of a cosmological constant (w ~ —1) if \d\og F/d\ogb\ <C 1. 
Given the FLRW metric (3.1), the matrix (2.5) reads 
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B = -^-(d<l>) T (I-v®v)d$ (3.5) 



a 2 

where I is the identity matrix, v <g) v is a matrix of elements v % v\ and we have defined 

(d®)° = di<$> a , v = -(d^)- 1 ® (3.6) 

where the dots denote derivatives with respect to conformal time r. Recalling that the total time 
derivative of $ vanishes, i.e. 

d§ a d$ a d$ a . , . 

-I- = + ^ v = 3 - 7 
dr dr dxi v ' 

one recognizes v in (3.6) as the (Lagrangian) velocity 

v i = x i (3.8) 

In the language of fluid dynamics, the equation (3.7) means that the material (or convective) 
derivative of $ is zero, which is nothing else than the statement that the label of a fluid element 
does not change along its trajectory, accordingly with the interpretation of $ as the comoving 
coordinates. This is equivalent to the first of the conditions (2.1). In fact, from the definition of 
the four- velocity, = dx^/drj, we get, consistently with equation (3.7), 

n dr 1 ,- dx l dx l dr 1 . 

u° = - = - 1 u* = — = —— = - 1 v t 3.9 
drj a drj dr dr] a 
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where 7 = (1 — v 2 ) 1 . 

The (square root of the) determinant and the inverse of (3.5) are then given by 

b = -idet(0$)Vl - v 2 (3.10) 

(X 

-1 Mr 1 1 *>> ra*T\-i 



B- L = a 2 d$- L I + -v <8> v (<9<T ) _1 . (3.11) 

\ 1 — v 2 J 

The equilibrium solution of the fluid in the FLRW background is given by = x l and therefore 

B lj = ^5 ij b=\ v l = 0. (3.12) 
a 2 a 6 

Moreover, the conserved currents can be explicitly expressed in a simple form 

J 5 )= ( p + pKW __ (3.13) 

^^ + ^(^[5} + —^ , (3.14) 

which is suitable for a Lagrangian formulation (see appendix B) of the conserved charges 

/f d 3 ® Fir 1 
fxa'jfan^ J - T (p + P)e a m wa u t (3.15) 

where n M is the unit normal vector to constant time hypersurfaces. As expected, the background 
carries, apart from the entropy, no charges: n^J^ = 0. The ground state configuration is thus 
neutral but supports charge excitations. 

Notice that among the various solutions of the transversality constraint (2.23) there exists a 
class [30], £ a (<E>) = a c e a b c db5 3 (§ — $) with <3> and a c arbitrary constants, such that the integral in 
(3.15) gives rise to a time independent exact 2-form Q (living in the internal three-dimensional 
<3>-space): 

n = dv v = 1 Ui^-d^ a q = q = o. (3.16) 

\ J o® a 

The conserved charges associated to this class of solutions of (2.23) are nothing but the coordinates 
of the dual 1-form: *Q = Q a d& a . In turn, the circulation V defined as the flux of Q over an 
arbitrary surface S in the internal $-space 

V= [n = I V a d<$> a = I (^-^] Ul dx l (3.17) 

Js Jds Jr=x(r,ds) \ ° J 

is conserved on-shell 

V = 0. (3.18) 

This result represents the generalization to FLRW metrics of the standard non relativistic Kelvin's 
circulation theorem. Notice that the factor {p + p)/b reduces to mn/s in the non relativistic limit 
for a system of particles of equal mass m and number density n, where s = b is the entropy density 
(as discussed in the previous section). 
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4 Perturbations 



In this section we present the basic set-up for the study of fluid perturbations within this formalism. 
We will now focus on fluid perturbations alone and include metric fluctuations in Section 5. For 
simplicity, we assume that there is a single fluid, since the generalization to an arbitrary number 
of non-interacting components is straightforward. The starting point is the background state 4 , 
$ 1 (t, x) — x\ to which we add small independent fluctuations n a in each space direction 

^ a = x a + 7r a . (4.1) 

The natural expansion parameters are tt and dn. Physically, this means that the actual fluid 
cannot be described only using the background comoving coordinates because there are small 
inhomogeneities at each point which effectively break the symmetry properties of the system. 
Recalling the introduction in Section 2, this assignment of coordinates corresponds to switching 
from the isomorphism between the system of coordinates O and the unperturbed background fluid 
J 7 , to another one in which the target space is now the actual imperfect fluid. 
We have 

<9$ = I + dn $ = tt- det 9$ = 1 + dirf - ^-drfdji? + ^(drf) 2 + ... (4.2) 
and expanding the determinant of the matrix B up to second order in n we get 

b = b(l + dtf - h 2 - l -d^ k d k ^ + \d^d^ + (4.3) 
where b = a~ 3 . Therefore, after eliminating total derivatives, the action for the perturbations is 

~g (4 2) + • • •) (4.4) 
with 

C^ = -^tf-&WY) (4.5) 

being the quadratic part of the Lagrangian for perturbations, where we have defined the speed of 
sound: 




c\ = . (4.6) 



4 There actually exist infinitely many background states = Aa;* defined by the constant (in space and time) 
parameter A which fixes the background value of b — A 3 /a 3 . Instead of keeping track of all the various factors of A 
by expanding around one of these vacua, e.g. by using <f>* = \(x l + tt^O, we simply work with A = 1 and recover 
the general expressions valid also for A =/= 1 by expressing the physical results in terms of the corresponding generic 
background values b, p = —F(b), Ft — Fb(b), etc. For example, two fluids described by the same function F but 
different backgrounds 6's, depending on the values of the A's, will support perturbations with different speeds of 
sound which, nevertheless, are still expressed by (4.6) that is just evaluated into different background values. 
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In the case of dust c s = 0, whereas for radiation the speed of sound reaches, as expected, the 
value c 2 s = 1/3. The no-ghost condition is simply bF b < 0. Recalling (2.14), this inequality is 
equivalent to p + p > 0. Imposing also that (? s > gives in turn F bb < 0. Notice that at this 
level, neglecting metric fluctuations, the dynamics of the perturbations (as well as that of the 
background) is entirely given in terms of F, F b and F bb . Higher order derivatives enter only in the 
self-interaction terms. 

Let us also point out that the contribution of Cfi to the action, the linear piece of the La- 
grangian £ n , vanishes because it gives rise to a total derivative that we can omit in (4.4), since 
the background solution depends just on time by construction. 

It is straightforward to link the four-velocity of the fluid to the 7r fields that we have introduced 
in (4.1). This can be done at any desired order in n by solving iteratively the equations (2.1), 
which are sufficient to determine the four- velocity up to a sign. At second order in n (and choosing 
u° to be positive) the result is 

u° = -(l+ -ir 2 ] u l = - (-tt* + n k d k 7r l ) . (4.7) 
a \ 2 J a v ' 

Notice that v, as defined in (3.6), can be expanded as 

v = —(I — dn)§ + . . . (4.8) 

which is consistent with the expression v l = au l , by comparison with (4.7). From the last 
expression we recover (4.7) by neglecting 7, which is irrelevant up to third order in n. 

In the following, we shall discuss the transverse and longitudinal modes of the Lagrangian (4.5) 
by decomposing the coordinate perturbation: 

TT = 7T L + 7T ± V • 7T ± = V X 7T L = (4.9) 

such that the second order Lagrangian for n can be written as 

CW = l(p + P)(nl + nl-cl(V.n L f) . (4.10) 

The generalization to iV fluids with different density functions F a (a — 1, . . . , N) is simple. 
Each fluid will be described by internal coordinates $^ which can be perturbed independently 
using different 7c l a . If we assume that the fluids do not interact directly among themselves, which 
means that the Lagrangian does not contain terms mixing different sets of $ a coordinates, the 
generalization of (4.10) is just 

1 N 

A/ 2) = g PT + W ^ Q - K± + *aL " <L& ' (4.11) 

a=l 

where we define the relative background energy density of each fluid as the following time function: 
Q a = pa/pr- As we will see later, when we introduce metric perturbations, the specific case in 
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which the n aL are the same for two or more fluids is of particular interest because it describes 
adiabatic perturbations between those species. We end this section writing the equations of motion 
for transverse 5 

^{a\p a +p a )7T aA _) =0 (4.12) 

and longitudinal modes 

*aL + (1 - 3 cl A )UnU - clAiy ■ v aL ) = i = {1, 2, 3} a — l,...,N. (4.13) 

In this equation c 2 aA represents the standard adiabatic speed of sound of the fluid with label a, 
which for this type of fluids is simply the ratio between the time variation of the background 
pressure and energy density: 

Pa = c 2 aA p a . (4.14) 

In general, the adiabatic speed of sound is defined in hydrodynamics as the derivative of the 
pressure with respect to the density. In our case: 



2 _ dPc 
ZaA ~ d Po 



(4.15) 



b=b 



This definition reduces to (4.14) for a fluid given by (2.7). Using the equations (2.13) and (2.14), 
it is straightforward to check that, in this case, the adiabatic speed of sound c 2 aA actually coincides 
with the speed of sound c 2 as that we introduced in (4.6): 

cl A = c 2 as (4.16) 

for any fluid of the form (2.7). This shows that the only possible speed of sound of a perfect fluid 
defined by the Lagrangian (2.7) is adiabatic. As we will see later, the longitudinal modes, which 
(in absence of metric perturbations) are given by the evolution equation (4.13), are intrinsically 
related to the adiabatic modes of the fluid. 

If we use the expression (4.3) to expand at first order the energy density and the pressure of 
a fluid defined by (2.7), we find that the corresponding perturbations are 

b F 

SmP = (p + P) V • 7T L S M p = (p + P) —ft~ V • 7T L (4.17) 

fb 

in agreement with (4.16). In these formulas we have introduced the subscript m to indicate that 
these are the perturbations coming exclusively from the variation of the matter part of the action 
(2.6). In the next section we include metric fluctuations to achieve a complete description of the 
perturbations. The equations (4.17) will turn out to be modified by a correction coming from the 
metric inhomogeneities; see equation (5.16). 



5 The equation (4.12) corresponds to the equation of vorticity conservation, see equation (3.16). In fact, at this 
order in n, we have *fi = —V x (a 4 (p + p)7T±). 
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5 Including metric perturbations 



Assuming that the theory of gravity is General Relativity 6 (GR) the full action takes the form 

S = SeH + Srn (5.1) 

where the first piece is the usual Einstein-Hilbert action 

1 



'EH 



16nG 



(Pxy/=gR (5.2) 



and the second one is the matter part of S, which for each perfect fluid of the type we have been 
studying is given by (2.6). 

Let us now focus on S m , with the aim of finding the terms that constitute the direct interaction 
between metric and matter perturbations. A straightforward way of doing this at any desired order 
is to expand the action S m using a functional series. We expand first with respect to the metric 
variations taking an arbitrary matter configuration 

S m [(p,g^ u + 5g^ u ] =S m [(p,g^ u ] + \ J dtxy/^gT^Sg^x) + ... (5.3) 

where ip globally represents all matter fields (e.g. dark matter, dark energy, etc) and T^ v is 
the gravitational energy-momentum tensor. Expanding now the equation (5.3) around a matter 
background we end up with the following action for the matter-gravity coupling at linear order in 
metric perturbations (but all orders in the matter fields) 



S n 



D I J d 4 x^5 M T^5g, v (x) (5.4) 



where SmT^ is the variation of T^ v induced by the matter perturbations. The computation of 
the matter-metric mixing terms is then very simple provided that we know the form of the energy- 
momentum tensor. In particular, for a perfect fluid (2.7), this is just given by (2.11). We will now 
see in detail how this works in the conformal Newtonian gauge. 

In the following part of this section we focus only on scalar perturbations that mix with 
longitudinal modes. Vector and tensor metric perturbations will be discussed in Section 5.1. In 
the conformal Newtonian gauge the perturbed FLRW metric is then diagonal: 

ds 2 = a 2 (-(1 + 2i))dr 2 + (1 - 2<f>)6 ij dx i dx i ) (5.5) 

If the universe did not contain any imperfect fluids at all, we would have ip = <j) at linear order 
in the equations of motion. However, we want to include the more general possibility that some 
fluids with anisotropic stress could also be present and therefore we will treat these two potentials 
as distinct variables. This, for instance, happens at very early times when the anisotropic stress 
of neutrinos cannot be neglected. 



3 The formalism can also be applied to theories of modified gravity. 
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Let us consider the energy momentum tensor of some species in their rest frame. Using that 
u 2 = — 1 , we find the following expression for the (lowest order in metric fluctuations) coupling 
between matter and metric perturbations (of scalar type) 

^^SuT^dg^ = -a\^5 M p + 30 <W) (5.6) 

where the fluctuations in the energy density and the pressure include all orders. Notice that for 
a perfect fluid defined by (2.7) this coupling gives the following contribution to the action for 
fluctuations at first order in the metric perturbations 

S m D- J d A xa\i) + 3c 2 s <P)p5m (5.7) 

where 8m = &mp/p denotes the intrinsic (not metric) relative energy density perturbation at all 
orders and c 2 s is given by (4.6). For linear matter perturbations of such a fluid we can use (4.17) 
and therefore the mixing between matter and metric perturbations is 

S m D S%L g = - J d 4 xa\p + p)(ij + 3^0) d t nl (5.8) 

which shows that metric fluctuations and transverse modes do not couple at this order. 

As a check of these results, we perform an explicit independent calculation by writing the 
square root b and y/—g in the perturbed metric (5.5). The matrix B has coefficients 

a 2 B ij = -(1 + 24))- 1 && + (1 - 20)- 1 <9 fc $ i <9 fe $ i (5.9) 

and the square root of its determinant is 



where we recall that the definition of the coordinate velocity is (3.6) and det(<9$) is still given by 
(4.2). In addition, the metric determinant is 

v^ = aV(l + 2^)(l-20)3 (5.11) 

so that the mixing term at linear order in the metric perturbations is 

y/^gC m D a 4 (-0 — 30) 5F + 3a 4 5(bF b ) (5.12) 

which, using (2.13) and (2.14), is in agreement with the formulas above. 

We can easily compute all (metric and matter) second order terms that come from the matter 
action (2.6) by simply doing a Taylor series expansion. The result is 

S<? =\jd , xa t [l+ w) p [i 2 - 4(d,7,' L f - 2«' + 3cj0) ay L ] 

+ i/dW,^ + 30. -3(1 + »)^ (5,3) 
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The first line in this expression corresponds to the mixing term (5.8) that we have just computed 
plus the piece that comes from (4.10), which is the purely matter part of the action that we already 
obtained by neglecting metric perturbations. The second line of (5.13) is the contribution of the 
matter action (2.6) to the action of the metric fluctuations. In order to get the full action for 
perturbations at second order we have to complete the metric part by perturbing the Einstein- 
Hilbert action (5.2). This can be found, for any gauge, in [48]. 

The generalization of (5.13) to iV fluids of this kind, interacting only via gravity, is straightfor- 
ward. We just have to sum the individual actions of the different components. Then, we can easily 
write down the equations of motion for the transverse and longitudinal modes of each component. 
The transverse modes are unaffected by the scalar metric perturbations at this order and their 
equations of motion are still given by (4.12): 

^ (a 4 (p a +pa)ir a ±) =0 a = l,...,N. (5.14) 

On the other hand, the equation (4.13) for the longitudinal modes gains a contribution from the 
metric perturbation, becoming: 

TT aL + H(l - 3c 2 Jn aL - c 2 Sa V(V • n aL ) - Vip - 3c 2 Sa V0 = a = l,...,N. (5.15) 

The expressions in (4.17) for the density and pressure perturbations get modified when metric 
fluctuations are included. Concretely, they are replaced by 7 

5 pa = {pa + Pa) ( V • 7T a + 30) 5p a = (? as 5p a (5. 16) 

with the speed of sound squared c 2 as defined exactly as before, in (4.6). Notice that both the density 
and pressure perturbations gain a term that depends on the metric potential <fi. In particular, for 
the total (matter plus metric) density perturbation of each fluid we write 

Spa = 5 M Pa + 3(p a + p a )(j) (5.17) 

where the matter part 8m Pa is given by the first expression of (4.17). 

The four-velocities of the fluids also change with respect to (4.7) due to the effect of the metric 
perturbations. The equations (2.1) are both still valid and we can again find the four- velocity 
solving them iteratively. Alternatively, proceeding in the same way that we have used to obtain 
(3.9) we get 

m° = ^7 ui = \^ vl 7" 2 = (1 + 2^) - (1 - 20)w 2 (5.18) 
where the coordinate velocity is still given by 

v' = = -r + ^0^ + ... (5.19) 
dr 



7 This expression for the energy density reminds of the Zcl'dovich approximation, Sp ~ V • n, in Newtonian 
Lagrangian perturbation theory (see for example [49]), with n playing the role of the displacement field. The extra 
metric perturbation term (f> in (5.16) comes from our relativistic formulation. 
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which is the same that we would derive from (4.7) (where metric perturbations are set to zero). 
This happens because v l is the solution of (3.7), which is the same with or without metric pertur- 
bations. Expanding 7 at second order in perturbations we get 




1 - rj) + -^j 2 + -v 2 + ... 



U* = -(i + 



(5.20) 



Let us point out that the expression (5.20) for u % differs (at second order) from the one found 
in the review [50] and other related papers in the literature. 

In the next section, we will use the equation (5.15) to understand adiabatic modes. In order to 
do so, we will need to replace and ip by their expressions in terms of pure matter perturbations. 
The difference between the two metric potentials is given, at first order, by the total anisotropic 
stress of the fluid system. If all its components are perfect, the metric potentials are equal at first 
order. Any difference between the metric potentials will then be due to imperfect components. 
Another useful piece of information is the (general relativistic) Poisson equation, which allows to 
express the Laplacian squared of as a function of the density and velocity perturbations. This 
equation can be derived by combining the — and — i Einstein equations (see for instance [14]) 
or from the variation of the second order perturbed action (5.1) with respect to a metric degree 
of freedom that is later set to zero in the conformal Newtonian gauge [48] . The Poisson equation 
reads: 



where we sum over all fluid species and 6 a is the standard notation for the divergence of the 
velocity perturbation v a [51] , as given in (5.19): 



Since for a perfect fluid, 5 a is given by (5.16), it turns out that the right-hand side of (5.21) 
contains a contribution that depends on the spatial metric perturbation 0. Hence the expression 
(5.21) becomes a partial differential equation for with second and fourth order derivatives. We 
will use such a form of the Poisson equation in Section 6 when discussing adiabatic perturbations, 
see equation (6.16). 

Before moving to the description of adiabatic modes, it is interesting and useful to see how 
the equations of motion and (5.14) and (5.15) are related to the standard continuity and Euler 
equations for generic fluids that have the following general form in the conformal Newtonian 
gauge [51]: 




(5.21) 



a=l 



9 a = V • v a = -V ■n a + ... 



(5.22) 




(5.23) 



e, 



(l-3w a )W, 



1 



a 



V 2 5 a - V V + VV a . 



(5.24) 



a 



a 



1 + W a 



l + w a 5p a 
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Notice that for these fluids, i.e. described by (2.7), the Euler equation (5.24) with a a = is 
precisely the divergence of (5.15), while the continuity equation (5.23) is just an identity. This 
follows from the fact that the velocity perturbations 6 a in (5.22) are the time derivative of the 
longitudinal modes V • n a rather than independent variables. 

5.1 Vector and tensor metric perturbations 

We have derived the equations of motion (5.14) and (5.15) for an ensemble of N non-interacting 
fluids, each with its own 7r a field, assuming only scalar metric perturbations. However, since n a ± 
can be written for each fluid as the curl of a vector potential, we can expect its dynamics to be 
affected by metric vector perturbations. For the longitudinal component n aL this cannot occur 
because it has zero curl and hence it can be expressed as the gradient of a single scalar degree of 
freedom, which then couples only to <fi and ip as we have already seen. 

The perturbed FLRW metric in the Poisson gauge [52,53] generalizes the conformal Newtonian 
gauge (5.5) to include vector and tensor degrees of freedom: 

ds 2 = a 2 (-(1 + 2i))dT 2 + 2i/ drdx i + [(1 - 20)<% + Xij] dx^x 4 ) (5.25) 

where v and x are respectively pure vector and tensor perturbations that satisfy 

d k v k = d kXkj = Xjj = 0. (5.26) 

Applying (5.4) to the full set of metric perturbations (transverse) v vector and (transverse and 
traceless) x tensor metric degrees of freedom we find 

^y^5 M T^5g, u =a A [-ty 5 M p + 30 5 M p) + a(p + <W] (5.27) 
= - a 4 (p + p) [(ip + 3c 2 s (j)) diiri + tt^z/] + (total derivative) (5.28) 

where tensor perturbations do not appear because of the tracelessness condition. This confirms the 
known result that tensor metric perturbations do not couple, at the lowest order in perturbation 
theory, to perfect fluids. On the other hand, the transverse modes (vortices) mix with the vector 
metric perturbations. The equations (5.14) get now replaced by a conservation equation of the 
linear combination (z/ — txt) 

^{a 4 (Pa+p a ){v t -K±))=V a = l,...,N (5.29) 

where, as before, the index a lists the different perfect fluids from 1 to N. This equation is nothing 
but the conservation of the charges associated with the vorticity, see (C.17). It corresponds also 
to the separate conservation of the three-momentum for each fluid [50,54]: 

5qi + mbqi = a — 1, ... ,N (5.30) 
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where 



Sqi = (Pa + PaW ~ r a± ) a = l,...,N (5.31) 

are the three- momenta of the fluids. Notice that, if the fluids were not perfect, the right-hand 
side of (5.30) (or equivalently (5.29)) would have an anisotropic stress source term. The equation 
(5.30) tells us that in absence of such a term, the three-momentum decays due to the Hubble 
expansion. 

The equation for the evolution of the vector metric perturbations is [54] 

V V = m 2 ]T(1 + w a )(S - r a± ) (5.32) 

a 

consistent with the redshift of the three-momenta at large scales. 

For completeness we include the equation for tensor modes, which in absence of any source of 
anisotropic stress reads: 

V 2 \„ = xa + WlXij ■ (5.33) 

As it is well known, these modes do not couple to any fluid degree of freedom if there is no 
anisotropic stress. 



6 Adiabatic modes 

Adiabatic perturbations correspond to modes associated with equal time shifts. They are therefore 
constructed by perturbing any homogeneous (intensive) fluid variable Ij in the following way 

Z?(t) ->■ Xj(r, x) = ij(r + 7r(r, x)) = Ij(r) + 7r(r, x)Ij + . . . (6.1) 

using the same 7r for all Ij , where the subscript j lists any intensive variable pertaining to the 
fluid. Therefore, we say that two intensive fluid variables I\ and X 2 are adiabatically related if 
their perturbations can be constructed from the same time shift n. If we focus on the energy 
density p and the pressure p of any fluid, at first order in n this implies 

c 2 A = c 2 s . (6.2) 

Moreover, for several fluids, we have that adiabatic density modes satisfy 

^ = ^ (6.3) 

Pa Pfi 

for any pair of fluids a and (3 interacting only through gravity. Describing the perturbations as in 
(4.1), we get that the adiabatic mode is nothing but a common longitudinal degree of freedom 

7r = ^ = ^ = --i-(V.7r + 30) Va. (6.4) 
p a p a SH 
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In other words, the condition of adiabatic modes translates into "flavour" (or species) independence 
of the longitudinal modes 

TTaL = 7T/3L • (6.5) 

In the next subsection we show how this result extends to imperfect fluids that have vanishing 
non-adiabatic pressure perturbations. 

For convenience, we use the curvature perturbation on uniform density hypersurfaces, which 
is defined to be [55] 

(^-t-H 5 ^. (6.6) 

Pa 

Remarkably, in our case this is just one third of the divergence of n a 

Ca = ^V-7T a (6.7) 

and taking its derivative with respect to conformal time, we get the following relation with the 
divergence of the velocity (at first order in ir): 

6 a + 3( a = 0. (6.8) 

Defining the entropy perturbation between two species in the usual way [56] 

S a p = 3(C* - Cs) (6-9) 

we see that for the kind of fluids that we consider, the entropy perturbation among two species is 
zero if the coordinate perturbations 7r Q (a — 1, 2) of the two fluids differ by at most a divergenceless 
three-vector. This is exactly the condition (6.5) for adiabatic modes that we have found before. 



6.1 Adiabatic modes for imperfect fluids 

Let us recall that for more general (imperfect) fluids 8 , entropy modes may come as well from an 
intrinsic non-adiabatic pressure perturbation. The total pressure perturbation of any species can 
be decomposed as 

5P a = 8P a{ nad) + c 2 aA 5p a (6.10) 

where the second term is the product of the density perturbation and the usual adiabatic speed 
of sound. For perfect fluids, 5P a ( na d) (the non-adiabatic part of the pressure perturbation) is zero. 
However, when a fluid has other internal degrees of freedom different from the $ coordinates 9 , 
an intrinsic non-adiabatic pressure will typically arise. Nevertheless, there exists an interesting 

8 For a study of cosmological perturbations of imperfect fluids, in particular in connection to scalar fields, see [57] . 
9 See Appendix A for the general form of the energy-momentum tensor in this case. 
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class of imperfect fluids that have vanishing non-adiabatic pressure perturbations. They fail to be 
perfect only because they have anisotropic stress a 

5>« + p a ) w a = J2 - \ E v2t " » ( 6 - n ) 

a i,j k 

that enters in the Einstein equations for the scalar potentials 

V 4 ^ - 0) = \u 2 J2 0.(1 + »a) VV Q . (6.12) 

At early times, neutrinos fall into this class of fluids: a v ^ while 5P u ( na( i) = 0. 

As we are going to see next, in an adiabatic mode, the expressions that occur for the density and 
velocity perturbations of perfect fluids are also valid at first order in 7r for imperfect fluids of that 
kind. From the equation (3.3) we know that S a /(l+w a ) is a species independent ratio for all fluids 
involved in an adiabatic mode. Since the definition (6.6) can be written as S a / (l + w a ) = 3(( a + (p) 
we have that 

Ca = C/3 (6-13) 

for any two fluids in an adiabatic mode. Moreover, if at least one perfect fluid is present, we have, 
thanks to equation (5.16), 3( = V ■ ir. For imperfect fluids with vanishing non-adiabatic pressure: 

or \ 1 + w a J 1 + w a l + w a 

where we have used the continuity equation (5.23) in the last equality. We thus find that in this 
case the velocity perturbations are also flavor independent 

SPa(nad) = 9 a = Of) = S£ (6.15) 

for any two fluids sharing an adiabatic mode. Again, if at least one perfect fluid is present, then 
9 is proportional to the time derivative of the longitudinal mode: 6 = —V • 7r. Effectively, in an 
adiabatic mode, the energy density and the velocity of an imperfect fluid with adiabatic speed of 
sound are like those of a perfect fluid. 

6.2 Superhorizon perturbations 

We now describe adiabatic modes in further depth using the variable ( introduced in (6.7). Since 
adiabatic modes are particularly important for the study of initial conditions in the early universe 
and at those times the anisotropic stress of neutrinos cannot be neglected, we are led to consider 
the generic situation in which the metric potentials <fi and ip can be different from each other, even 
at first order in the equations. Armed with the knowledge of the previous sections, we can write 
the Poisson equation (5.21) in the following way: 

(V 2 - M 2 ) V 2 = M 2 (V 2 C + 3?*C) (6- 16) 
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where we introduce the "effective mass" (squared) 

M 2 = 9 -n 2 J2^(l + w a ) (6.17) 

a 

in which the sum extends over all (perfect and imperfect) fluids that are present in the system. 

According to the Euler equation (5.24), for fluids with anisotropic stress but zero non-adiabatic 
pressure, the divergence of the equation (5.15) must be replaced by 

4 + (1 - K)UCa - <V 2 C - \v 2 ^ - <V 2 + \v 2 a a = (6.18) 

which, together with equations (6.16) and (6.12), form a close system for the perturbations. If 
photons (which have c 2 s = 1/3 and a 7 = 0) are present, recalling that for adiabatic modes the 
equation (6.13) holds, we obtain that the equation (6.18) implies 

C-^V 2 (C + ^ + 0)=O (6.19) 



3c 2 . 



m + \v 2 (c + 0) 



\v 2 a a =0 . (6.20) 



These two equations are actually valid only approximately. The same happens for the condition 
(6.13), which strictly holds just in the limit of large scales, V/"H — > 0. At early times, when 
the universe is radiation dominated and the anisotropic stress of neutrinos cannot be neglected, 
the equation (6.20) implies that V 2 a u = 0. Using the equation (6.12), one would then naively 
conclude that for an adiabatic mode V 4 (^ — (f)) = . These two last conclusions have to be 
understood as approximately valid expressions in the limit V/"H — » 0. Otherwise they would 
imply trivial solutions for o v and the difference of the metric potentials [51]. In the limit of very 
large wavelengths, the Poisson equation (6.16) simplifies to 

_V 2 =3U( + V 2 C • (6.21) 

If we plug (6.21) into (6.20) (which assumes adiabaticity), we obtain V 2 a a = 0, without having 
made any assumption about c 2 a . This means that the approximation V 2 a u = that we have 
found for adiabatic modes is valid at very large scales (and early times, with the universe being 
radiation dominated). Therefore, as we have anticipated, we conclude that adiabatic modes can 
only be defined at very large scales and, strictly speaking, there are no exact adiabatic modes. 

Notice that if the perturbations are regular for very large wavelengths, the Poisson equation 
helps us to write (6.18) in the following way 

C + X - V 2 C + 2H( = *-U 2 + w a )a a (6.22) 

a 

which is a simple equation for ( where no metric perturbations appear explicitly. If the condition 
of regularity at large distances wouldn't apply, the whole equation (6.22) takes an overall extra 
V 2 operator on both sides. 
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To conclude this section, let us recall and show that, as it is expected, the curvature per- 
turbation on uniform density hypersurfaces ( coincides with (minus) the comoving curvature 
perturbation 1Z in the limit of very large scales. The definition of 1Z is (see e.g. [47]): 

K = <l>-Uv (6.23) 

where v is a scalar potential that gives the longitudinal part of the fluid coordinate velocity, i.e. 
n L = —Vv . Taking the Laplacian of (6.23) and using (6.7) we get 

V 2 ft = V 2 + 3HC (6.24) 

For perturbations of very large wavelength, we can use (6.21) to eliminate the time derivative of 
( in (6.24), obtaining 

Tl + ( = (6.25) 
assuming again regularity of the fluctuations at very large distances. 



7 Conclusions and Outlook 

We have shown how to obtain the fundamental equations for the cosmological evolution of perfect 
fluids from symmetries and action principles. Although we have chosen to focus our analysis 
on the case of a FLRW metric (and its perturbations) in General Relativity and non-interacting 
fluids, the method is of ample generality and can be easily applied to other geometries, theories of 
modified gravity and can be extended to describe systems which do not interact only via gravity. 
The formalism presented allows a straightforward description of the dynamics of the expansion 
both at the background and perturbation level. 

The method is based on an effective theory which is an expansion in spacetime derivatives 
of three scalar fields. The condition of invariance under spacetime diffeomorphisms and inter- 
nal volume preserving spatial diffeomorphisms is sufficient to characterize completely the formal 
structure of the theory. At the lowest order in derivatives, the action for perfect fluids is given 
by just an arbitrary function F of a single operator that takes the form of a determinant be- 
cause of symmetry requirements. Identifying the various thermodynamical quantities, it turns out 
that F is minus the energy density and is a function of the entropy density. Its first and second 
derivatives are enough to characterize the pressure, the temperature and the speed of sound of the 
fluid. Higher order derivatives enter in the self-interactions of the perturbations. Different types 
of fluids can thus be described by adequate choices of the functional dependence of the energy 
density on the entropy density. Standard cold dark matter and radiation, just to mention two 
simple important examples, are easily described by specific power functions. 

Considering the high precision of ongoing (e.g. [58]) and upcoming (e.g [59]) observations, the 
need for an accurate theoretical understanding of the Universe's content is more important than 
ever. In order to interpret the data appropriately, a framework that describes the cosmological 
dynamics of the matter (energy-momentum tensor) in wide generality can be a leap forward. A 
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natural and widely popular approach has been that of using general fluids (perfect or not) [60] 
in a phenomenological and quite ad-hoc fashion to fit different types of data. As an important 
example, many of the studies that try to determine the properties of dark energy are based on 
the choice of a fluid with a (possibly time- varying) equation of state close to —1 and some model 
dependent assumptions for its perturbation properties. Conversely, the flexibility of the formalism 
presented here can be relevant for model comparison and parameter estimation since the physical 
properties of a perfect fluid are encoded in a single function that defines the action. Instead of the 
common phenomenological approach of fitting several parameters like relative densities, equations 
of state and speed of sounds of an ensemble of fluids, one could directly work with the functions 
F's that define the fluids. Therefore, we think that the framework presented in this paper can 
provide a first step towards a different useful approach in which the search is focused on functions 
F that fully describe the energy momentum tensor. 

The interaction terms between matter (at any order) and metric (at first order) perturbations 
can be easily extracted from the gravitational energy-momentum tensor. This is the case for 
scalar modes but also for vector and tensor ones, see equations (5.7) and (5.27). The knowledge 
of the function F and its first two derivatives is enough to describe perfect fluids at first order in 
perturbations around the background solutions. Higher derivatives of F will appear at the next 
orders in matter perturbations, which can be easily incorporated by using the simple structure of 
(5.7) and (5.27). 

The longitudinal (compressional) modes are found to be intimately related to the adiabatic 
modes at early times, ( — V • 7r/3, while the Goldstone boson 7r of the associated time shifts 
contains a piece that depends on the metric perturbations, 7r = —(( + 0)/"H. We have shown 
that the equations that govern the evolution of adiabatic modes (6.19)-(6.22) remain true even 
for non-perfect fluids which have vanishing non-adiabatic pressure (such as neutrinos at early 
times) provided that at least one fluid species, e.g. photons, can be described by a perfect fluid. 
Transverse modes (vortices) mix with the vector metric perturbations and evolve accordingly to 
the dilution of vorticity under Hubble expansion. 

There are several directions that are worth exploring in further depth. Among these, we find 
particularly interesting the possibility of adding more symmetry such that the system is more 
constrained and predictive. Conversely, demanding less restrictive symmetries would allow to 
draw more general conclusions. Below we discuss two such examples that are relevant for almost 
de Sitter expansions. 

More symmetry: scale invariance 

The condition of almost de Sitter expansion, either at early times (inflation) or at late times (dark 
energy), is w ~ —1, that is bF^/F <C 1, see equation (2.17). It corresponds to the condition that 
the action is weakly depending on b, suggesting the idea of adding extra (approximate) symmetries 
that could suppress F(6)-type of actions 10 . An example of such a symmetry is provided by internal 

10 However it is important to stress that if the would-be leading term F in the derivative expansion is suppressed, 
then next order terms in the derivative expansion should generically be taken into account. 
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scale transformations 

$ a ^e A $ a , b^e 3X b. (7.1) 

Notice that the choice F(b) = /tin b gives an action that is invariant only up to shifting the 
cosmological constant, F(b) — > F(b) + 3kA. The dynamics associated with such an action would 
thus be symmetric under (7.1) only if k <C 1 or gravity is switched-off. 

Other symmetries may also be invoked to sort the desired Lagrangian F(b) among all possible 
choices. An interesting case is e.g. Weyl symmetry 

9^(x) -+ e^g^x) $ a (x) -+ <5> a (x) (7.2) 

where u is an arbitrary function of the spacetime coordinates. $ has vanishing scaling dimension, 
in agreement with the fact that the identically conserved entropy current must have scaling di- 
mension 3. Scale transformations (and the full conformal symmetry) is a subgroup of Weyl x Dif- 
feomorphisms. Under Weyl transformations b — > e~ 3uJ b and thus the action is invariant when 
F(b) oc b A l z which corresponds to the equation of state p — |p of radiation. Such a fluid has 
indeed vanishing from the conservation of dilation current. Notice that the diagonal (internal 
and spacetime) scale transformation A = wA =constant 

g^(x) -+ e 2 "g^(x) $ a (x) -+ e" A $ a (x) (7.3) 

that gives b — > e~ 3 ^~ A ^ u b leaves the action invariant when A ^ 1 and F(b) is a homogeneous 
function of rank 4/[3(l — A)]. Radiation and matter are selected requiring invariance under (7.3) 
with A = and A = —1/3 respectively. 

Less symmetry: solids 

It is interesting to consider the possibility of reducing the symmetry of the system by demoting 
the volume-preserving internal diffeomorphisms (2.4) to be an approximate symmetry. This would 
correspond to a homogeneous and isotropic solid, i.e. a jelly, with very small transverse speed of 
sound c\ <C 1. For a jelly the dynamics is still invariant under translations (2.2) and rotations 
(2.3) and the action depends not only on b but also on the traces of B and B 2 

C = F{b, Tr B, Tr B 2 ) . (7.4) 

Assuming e.g. C = F(b) + eTrB one generates a 0(e) speed of sound for the transverse modes ir±. 
This could be useful within inflationary models when trying to quantize the perfect fluid since it 
allows to damp the wild quantum transverse fluctuations around the classical vacuum = x % [31]. 

Note added 

After completion of this manuscript, a related work treating the effective theory of solids for 
inflation appeared [61]. 
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Conventions and notation 

We use natural units so the speed of light is simply c = 1 and the reduced Planck constant is 
h = h/(2ir) = 1. Our choice of Lorentzian signature is mostly positive: (—,+,+,+). Tensor 
components in 4-dimensional spacetime are denoted by Greek indices running from to 3, with 
corresponding to the time coordinate. Latin indices, running from 1 to 3, are used for the purely 
spatial part of a 4-dimensional tensor or quantities intrinsically defined in M. 3 (such as T^..) or in 
the internal fluid space T (e.g. T a b c ...). In any case, Greek or Latin, repeated indices in a product 
will be summed over by convention, unless it is otherwise specified. The symbol V denotes the 
3-dimensional spatial gradient. Instead, 4-dimensional covariant derivatives are denoted with a 
semi-colon subscript. For example, the covariant derivative of a 4- vector is denoted by n M ;M . We 
use € ! j_ U p cr to define the 4-dimensional totally antisymmetric Levi-Civita symbol with €0123 = 1- 

The 

symbol e^ upa is totally antisymmetric and has e 0123 = — 1. The 3-dimensional totally antisymmetric 
Levi-Civita symbol is = e^jk- 

We distinguish a background quantity form its full inhomogeneous counterpart by denoting 
the first one with an overbar. For instance, p denotes a background (homogeneous) density while 
p denotes the corresponding inhomogeneous variable. 

Partial derivatives with respect to conformal time r are denoted with overdots. For example 
H = a/ a is the conformal Hubble parameter. 

A Relativistic hydrodynamics and the energy-momentum 
tensor 

We now interpret the fluid defined in Section 2, equation (2.7), in the context of standard general 
relativistic hydrodynamics. The scope of this section is understanding how restrictive are our 
assumptions of homogeneity, isotropy and invariance under volume preserving diffeomorphisms of 
the internal <3> coordinates. As it is well known, see e.g. [62], any energy-momentum tensor can 
be expressed in the following covariant way 

= p u^Uy +ph„ v + q u + u u q^ + 7i> (A. 1) 
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where is a choice of frame (a timelike four-vector) and h^ u = g^ u + u^u v is a projector on 
hypersurfaces orthogonal to . The energy density, pressure, heat flux and anisotropic stress are 
respectively defined as the following projections of the energy momentum tensor: 

p = T flu uV 3p = T flu h^ % = -hlT v y n^ = hffiT a p-h IIU h af jT a P. (A.2) 

By construction the heat flux is orthogonal to the four-vector with respect to which we decompose 
the energy momentum tensor: = 0. Assuming the weak energy condition 

T^U"U U > fW M = -1 U° > (A.3) 

(for all future pointing timelike four- vectors £/ M ), there is a unique timelike eigenvector u?™ of 
T M1 , that is of unit norm, u 2 E ^ = — 1. This four- vector defines the so called energy (or rest) frame 
of the fluid [63]. In this frame the heat flux is zero, q? E * = 0. For a perfect fluid (or one that is 
in thermodynamical equilibrium) the particle flux and the entropy flux are both parallel to 
which thus defines the unique choice of hydrodynamical four-velocity for the fluid 11 [39,41,62]. 
This can be expressed as 

u" = ^ (A.4) 
drj 

where dr\ = \J—ds 2 is the proper time. The energy momentum tensor of a perfect fluid (i.e. with 
TT/ii/ = 0) is then 

= (p + P) U(E)nU( E )u +pg^v (A.5) 

When dealing with perfect fluids we drop the subscript (e) to simplify the notation. From (A.5) 
we see that imperfect fluids in their own rest frame display anisotropic stress, which can originate 
from dissipation or non-gravitational interactions with other fluids [64]. 

Perfect fluids like (A.5) are barotropic if the energy density and pressure are functionally 
related, p = f(p). Non-barotropic perfect fluids still admits an equation of state but the energy 
density (or pressure) in the rest frame does not completely characterize the thermodynamical 
properties of the system. The perfect fluid defined by (2.7) is barotropic since both p and p are 
functions of the single variable b. 



B From Euler to Lagrange 

In this appendix we comment further on the relation between the Eulerian and the Lagrangian 
formulation of the dynamics of perfect fluids. Going from the proper time dr/ = \/—ds 2 to the 

n If the entropy or particle flows pointed in other directions, it would be possible to define other frames, following 
the world lines of these flows. 



24 



conformal time r, and from the comoving coordinates $ to the physical space coordinates x, by 
means of equation (3.10): 

^d 4 x= d -^dr ] (B.l) 

in the FLRW background metric. Actually, this result is valid for arbitrary metrics, as it is clear 
from (C.4) and (C.12) in Appendix C. The action (2.6) can thus be rewritten in a Lagrangian 
formulation 



J d A x^gC m = - J d 3 <5>d V ^, (B.2) 

where we have used equation (2.13), which tells us that the Lagrangian C m = F is equal to 
the comoving energy density p = T^u^u" multiplid by —1. Notice that d 3 Q/b is the invariant 
3-volume form on spatial sections at constant times. This means that 

&(*) " b(f) ' (B ' 3) 

where <3> a — > / a ($) is a 3-dimensional diffeomorphism. 

We see from equation (B.2) that the action for the fluid is just the continuum limit of a 
(iV-body) "mechanical" system made of point-like particles carrying no charges and being char- 
acterized only by their trajectories in spacetime 



S 



~^2 m i / dl li — > S m = - drj—p. (B.4) 



N 



Clearly, the continuum limit has the advantage, over the corresponding discrete (and relativistic) 
iV-body problem, that gravitational interactions and the matter backreaction are easily taken into 
account by promoting the metric to a dynamical variable and performing a local perturbation 
theory as presented in Section 5. Moreover, while iV-body problems deal with 3iV 1-dimensional 
degrees of freedom, the continuum limit describes the system in terms of just 3 degrees of freedom 
in a 4-dimensional manifold. The conservation of vorticity renders trivial the dynamics of 2 degrees 
of freedom. 



C Vorticity conservation 

In this appendix we use the ADM formalism [44] to provide the expressions (valid for any metric) 
of the conserved charges associated with the internal volume-preserving diffeomorphisms (2.4). 
In the ADM parameterization, the metric reads 

ds 2 = -N 2 dr 2 + {dx l + N'drfhijidx* + N 3 dr) (C.l) 

and its inverse has components 

a 00 - 1 a™-*' aV-W- NiNj (C2) 
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The determinant of the metric is given by ^/~# — Ny/deth. We define the three-dimensional 
vector 

e = + AH, (C.3) 

in terms of which the fluid four- velocity = dx 11 / drj (being r\ conformal time) is 

u° = i , u l = ^ (C.4) 

where 

= -Kd^r 1 ]^ = -— (c.5) 

and, as in Section 3, <9$ is the matrix of elements ((9$)/ = d^ a . The scalar product £ 2 is 
computed according to the induced spatial metric h, 

e = e^Kj . (c.6) 

This allows to write concise expressions for the matrix B defined by (2.5), its inverse B' 1 and 
the square root of its determinant (2.8): 

B = <9$ T [h- 1 - £ <g> f] 0$ , (C.7) 
B- 1 = (d®)- 1 [h + u®u] (d$ T ) _1 , (C.8) 



b = detm ^^t ' (C,9) 

The outer products that appear in the expressions (C.7) and (C.8) are [£ ® = and 
[u®u]ij = u,iUj with Ui = g^u^ = / y/l — £ 2 . The inverse spatial metric, that we denote h~ l , 
has elements h 1 ^ (such that h^hjk = S l k ). The spatial indexes of purely three-dimensional vectors 
are lowered by the spatial metric, = and iVj = h^N^ = g 0i . 
As discussed in Section 2, any function e a ($) gives rise to a current 

J£ } = -6F 6 e c ($)(S- 1 ) 0,, ^ d (CIO) 

which is conserved on the equations of motion (2.9) as long as d a e a = 0. The transversality 
constraint (2.23) ensures that $ a — > $ a + e a ($) is, locally, a volume-preserving diffeomorphism. 
The associated conserved charges are defined in the usual way 

Q {£ ) = - J (FxVd&hJfcnr Q {£) = 0, (C.ll) 

where n M = 1/A r (l, —N % ) is the normal unit vector to the constant time hypersurface S r that has 
h as the induced spatial metric. 
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Using (C.9) to pass to the Lagrangian formulation: 

b Vdet h 

the conserved charges can be expressed in a simple form 



Q(e) = ~ 



(P + P) a 9 ^ 



(C.12) 



(C.13) 



Solving the constraint (2.23) with £ a ($) = e abc d b f c (&) the charges take the following form after 
integrating by parts: 



f a ($)e abc d b 



(p + p) dx % 
b d$ 



-Ui 



(C.14) 



Let us now move on to the specific case of the vorticity currents, generalizing the results of 
Section 3. The vorticity charges are defined by picking a delta function in (C.14), f a = a a 5 3 (&— <E>). 
The resulting vorticity charges are given by the curl, in the three-dimensional <3>-space, of a 
(co)vector V 



Qa — £abcd b V c , 



(p + P) dx i • 

V a= ; T^-Ui, Q a = 0, 



(C.15) 



b d$ a 

where we have renamed $ — > $. By comparing these expressions with (C.14), we see that any 
conserved charge can be expressed as a linear functional of the vorticity charges Q a : 



Q U) = J rf 3 $/ a ($)Q a ($) . 



(C.16) 



Notice that the leading order term in perturbations for Q, around the metric (3.1) with the 
(background) field $ a = x a choice, gives 



Qa = ~ 



a 4 (p + p)e aij di tt j ± - 



+ 



(C.17) 



where the ellipsis contain all the higher order terms. This shows explicitly that in the Poisson 
gauge (5.25), the vorticity conservation is precisely equivalent to the dynamical equation (5.29) 
for the transverse modes. 

The common definition of vorticity circulation V can now be naturally extended to arbitrary 
spacetime metrics: it is the flux of the vorticity Q over a surface S (with boundary dS) in the 
internal $-space 



y = f *(Q a d$«) = I V a d$ a = I 

JS JdS JV 



(P + P) 



Uidx' 1 



(C.18) 



where the boundary OS is mapped into a closed loop, V = x(r,dS), in real space. Since Q is 
constant on-shell, the circulation V is conserved on the solutions of the equations of motion. This 
shows that the results of [30] for Minkowski spacetime carry over arbitrary metrics that admit a 
3+1 decomposition. 
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